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1 Materials and Methods 

1.1 Samples 

High-quality single crystals of Pbo.ssBii sSri 6Lao.4Cu06+5 (Pb-Bi2201) and 
Bi2_3;Pba;Sr2CaCu208+5 (Pb-Bi2212) near optimal doping were grown by the travelling-solvent 
floating zone method. The Pb doping suppresses the super-modulation in the BiO plane and min- 
imizes complications in the electronic structure due to photoelectron diffraction. The carrier con- 
centrations of the samples were carefully adjusted by a post-annealing procedure under flowing 
A^2 gas which varies the oxygen content, and were estimated from the Fermi surface volumes, e.g., 
for Pb-Bi2201 shown in Fig. 1 to be ~ 25.3 %, which is consistent with previous reports for near 
optimal doping (SI, S2, S3). X-ray and Laue diffraction showed no trace of impurity phases. For 
Pb-Bi2201, the onset temperature of the superconducting transition, Tc, determined by SQUID 
magnetometry, was 38 K, with a resistive transition width less than 3 K (Tc = 97 ± 2 K for Bi2212 
related to Fig. 4D); the pseudogap temperature, T* = 132 K ±8 K, is determined based on the 
observed closing of the antinodal pseudogap in ARPES (Fig. SdP), which coincides within exper- 
imental uncertainty with the onset of Kerr signal in PKE and that of transient reflectivity in TRR 
(Fig. 3). This T* value is also consistent with our resistivity measurement on the same crystals as 
well as other reports of different measurements on Bi2201 near optimal doping (S3, S4, S5, S6, S7). 
We found T* depends more sensitively on the post- annealing condition than Tc near optimal dop- 
ing but remains consistent among samples within the same growth and post- annealing batch. 

1.2 Measurements 

1.2.1 ARPES 

Angle-resolved photoemission spectroscopy (ARPES) measurements were performed at Beamline 
5-4 of the Stanford Synchrotron Radiation Lightsource (SSRL) with a SCIENTA R4000 electron 
analyzer. Preliminary photon energy and polarization dependent ARPES measurements were per- 
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formed at Beamline 10.0.1 of the Advanced Light Source (ALS). All presented data were taken 
using 22.7 eV photons mainly in the first Brillouin zones with total energy and angular (momen- 
tum) resolutions of ~ 10 meV and ~ 0.25° (~ 0.0096A~^), respectively. The temperatures were 
recorded closest to the sample surface position within an accuracy ± 2 K. The samples were 
cleaved and measured in an ultra high vacuum chamber with a base pressure of better than 3 x 10~^^ 
Torr which was maintained below 5 x 10~^^ Torr during the temperature cycling. Measurements 
on each sample were completed within 48 hours after cleaving. All data shown in the main text 
were obtained on the same sample and were reproduced on different samples from the same batch 
with sample cleaving at low (10 K) or high (~ 150 K) temperature. The possibility of sample 
aging in each experiment was excluded by monitoring the nodal dispersion without complications 
by the gap opening. A similar experimental routine has been previously demonstrated and detailed 
in the supplementary information of Ref. (S3). 

1.2.2 Polar Kerr effect 

Polar Kerr effect (PKE) measurements were performed using a zero-area-loop Sagnac interferom- 
eter at a wavelength of A = 1550 nm, spot size of ~ 3/im, and sensitivity of 0.1 /irad/\/ Hz at 250 
/iW of incident optical power (S8). In its used configuration, the apparatus was sensitive to only 
the polar Kerr effect, hence to any ferromagnetic component of the local magnetization, perpen- 
dicular to the plane of incidence of the light. The same apparatus was previously used to detect the 
effect of time-reversal symmetry breaking below Tc in Sr2Ru04 (S9) and in the vicinity of T* in 
YBaaCugOe+a: (YBCO) (SIO). 

For YBCO, measurements on crystals with various doping levels revealed a sharp phase tran- 
sition at a temperature T's(p) below which there is a non-zero Kerr rotation. Both the magnitude 
and doping dependence of Tg were found to be in close correspondence with T* which has been 
identified in other physical quantities. In particular, Tg is substantially larger than Tc in under- 
doped materials, but drops rapidly with increasing doping, so that it is smaller than Tc in a near 
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optimally-doped crystal and extrapolates to zero at a putative quantum critical point under the su- 
perconducting dome. The magnitude of the Kerr rotation in YBCO is smaller by ~ 4 to 5 orders 
of magnitude than that observed in other itinerant ferromagnetic oxides, suggesting that at most 
it measures a very small "ferromagnetic-like" component of a magnetic transition. However, the 
temperature dependence of the PKE near T^, together with complementary measurements using 
other probes, strongly suggest that the PKE tracks a secondary order parameter which is driven by 
another transition that is not necessarily magnetic. Corroborating evidence in that direction comes 
from neutron measurements on x= 0.45 (Sll) finding evidence for a nematic state in this composi- 
tion, and from muon- spin-rotation measurements indicating charge ordering transition (572) near 
optimal doping. Interpolating between these two ends which find T* that agrees with the PKE 
measurements, we speculate that the PKE tracks an electronic transition that may evolve smoothly 
from strong charge ordering to a weaker nematic phase. 

While the above discussion on YBCO relies on circumstantial evidence, the present results 
on Pb-Bi2201 may be the first unambiguous data that suggests that T* can be identified, which 
corresponds to a true symmetry-breaking transition and the primary order parameter that governs 
that transition is electronic. 

1.2.3 Time-resolved reflectivity 

Time-resolved reflectivity (TRR) measurements were performed using a 100 MHz mode locked 
Ti: Sapphire oscillator operating at a wavelength of 800 nm. In these measurements, femtosecond 
(fs) pulses of linearly -polarized light excite the sample, and time (t)-delayed probe pulses from the 
same laser measure the resulting change in reflectivity, AR. In the present context, the power of 
this technique lies in its ability to distinguish different phases by their response dynamics. Fig. 
S(3] shows TRR data taken on Pb-Bi2201 at 3 different pump powers with a focal spot diameter of 
roughly 100 fim. The grey and pink shaded regions highlight the pseudogap and superconducting 
responses, respectively. As discussed in the main text, the superconducting signal, seen in Fig. 
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S|3^ & C, is attributed to the breaking of Cooper pairs, and its temperature and pump-fluence de- 
pendence are consistent with previous work {S13). It is characterized by a positive signal which 
decays on the picosecond (ps) time scale. The pseudogap signal, on the other hand, is centered 
about zero delay and decays on the 100 fs timescale. The pseudogap signal is generally weaker 
than the superconducting one and becomes more pronounced with increasing pump power. How- 
ever, the average heating of the sample becomes significant at low temperatures for high pump 
powers, which can almost completely suppress the superconducting response. The magnitude of 
the negative TRR signal in the left inset of Fig. 3 is extracted from the 15 mW data (Fig. S|3}\) by 
taking the difference between the maximum for t < and the minimum near t = 0. Because of 
the laser heating at this power, we only show the data for T > 25 K, and estimate an error bar of 
±8K for the actual sample temperature. 

1.3 ARPES data analysis 

To remove the effect due to cutoff by the Fermi-Dirac (FD) function, raw ARPES spectra have 
been divided by a convolution of the FD function at the given temperature and a Gaussian with its 
width fixed at the energy resolution. Note that this division is an approximation for the complicated 
deconvolution that aims at removing the energy resolution effect. It works well for high tempera- 
tures when the overall width of FD function (~ AksT) is much larger than the energy resolution 
and becomes bad for low temperatures. This procedure allows us to recover the actual band disper- 
sions closest to Ep and trace them above Ep at high temperatures, when thermal population leads 
to appreciable spectral weight. It is also necessary for revealing that the EDC shoulder feature at 
low temperatures is not trivially produced by the Fermi cutoff. Because the uncertainty introduced 
by the division to the spectral line shape is limited to ~ 4 meV (= AksT) at 10 K around Ep, 
the discussed low-temperature spectral evolution across Tc for features (such as the EDC shoulder 
feature) located beyond this Ep vicinity remains robust. The faithfulness of this procedure can be 
reflected by the observations after the division that the EDC shoulder feature at low temperatures 
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does not show up everywhere but is confined in the antinodal region and that the nodal Fermi cross- 
ings are recovered at all temperatures. However, one should exercise caution when associating the 
gap function at low temperatures directly with the peak position of the FD-divided spectrum (Fig. 
$4]), particularly in the vicinity of the node where the gap is intrinsically small. 

The subsequent EDC analysis was performed on the FD-divided spectra: the EDC maximum 
below Ep is identified by taking the first derivative of the moderately-smoothed EDC (w.r.t. the 
energy axis); suspected additional features below Ep and local maxima above Ep are estimated 
based on the first derivative and aided with eyeball correction; the (non-)existence and energy po- 
sition of the EDC shoulder feature are determined mainly based on a spectral subtraction analysis 
(Fig. SO, which is in general consistence with the spectral division method using 40 K data (Fig. 
4C) and with eyeball estimates. 

2 Simulations 

2.1 Simple density- wave pseudogap order 

We describe modifications of the band structure due to (i-wave superconductivity (Fig. $7^ & 
F) and its coexistence with some long-range density wave order of wave vector q (Fig. S|7]C!, D, 
G & H) in a simple mean-field approach, modified from what we previously used in Ref. (S3). 
We confine the use of this approach only to low temperatures, where order parameter fluctuations 
are less severe and thus a simple mean-field description of the data might have a better chance of 
success. 

The mean-field coexistence Hamiltonian is given by: H = Y.k^kc\ck + Y.k,q^qi.A+q^k + 
h.c.) + Y.k ^k{c\-iC^-ki + h.c), where Vq and A^. are the density wave and superconducting or- 
der parameters, respectively; c|,(ca,.) is the creation (annihilation) operators for electrons at k. 

Ek = —2t{coskx + cos ky) — 4t' cos kx cos ky — 2t"{cos2kx + cos2ky) — At'" [cos 2kx cos ky + 
cosk:cCos2ky) - eo, where t,t',t",t"',eo = 0.22,-0.034315,0.035977,-0.0071637,-0.24327 
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eV, respectively, is the tight-binding bare band dispersion which is obtained by a global fit to the 
experimental EDC dispersions at 172 K (Fig. S|5]). 

We consider in the following only the density wave orders with generic gs that largely connect 
the antinodal portion (rather than the nodal one) of the Fermi surface, where the particle-hole 
asymmetry of spectra is the strongest. 

In case of bond-diagonal density wave order with a commensurate q = qaf = [tTjTt], the 
eigenstate {ipk) = Uk\k) +Uk+qAp\k + qAF) +U-k\-k) + 

U—k—qAF I — ^ — Qaf) for the Hamiltonian 

Aaf Dsc 



with the eigenenergy e'{k) can be obtained by solving the matrix Maf+sc 
where Aaf = 



Dsc —^AF 

Dsc = / is a 2 X 2 identity matrix. We expect this case to 



Sk V2 

generically represent density wave orders of q approximately along the bond-diagonal direction. 



such as (i-density wave, antiferromagnetic order with q = qAF and incommensurate spin stripes 
with a finite, small deviation of q away from qAF- 



In case of bond-direction density wave order, we consider the checkerboard case with qi 

^CB Dsc 
Dsc —AcB 

1600 X 1600 (identity) matrices. Acs is determined by the following, base: {k+mqi +nq2] m, n 



[0.157r,0] and q^ = [O.O.lSvr], Mcb+sc 



where Dsc = ^I, ^cb{I) are all 



0, ...,A^- 1; = 200/gcd{2* 100,0.15* 100) = 40}; matrix element: AcBim,n;m' ,n') = Vi, 
only if \m — m'\ + \n — n'\ = 1 or {\m — m'\ + \n — n' — N + l\){\m — m' — N + 1\ + \n — n'\) = 
[with small high-order interaction neglected {S14)], and otherwise. Eigenstate is thus \'ipk) = 
J2m,n Uk+mqi+nq2\k+mqi+nq2) +u_k-mgi-nq2 \-k-mqi -ngg)- Superposition of two orthogonal 
stripe orders with qi and q2 gives overall similar results as the checkerboard case. 

Several key aspects of our experimental observations (Fig. $7K & E) are qualitatively repro- 
duced in both these cases. In the nodal region, a single dominant branch opens up a rf-wave gap 
along the underlying Fermi surface (Fig. S|7B-H and the insets). In contrast, multiple branches 
of comparable spectral weight are seen in the antinodal region, with one branch exhibiting little 
dispersion in two dimensions and another dispersive branch which shows back-bending at kc away 
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from kp. 

2.2 Pair-density-wave pseudogap order 

Other density-wave alternatives for the broken-symmetry pseudogap order have not been (but 
should be) explored. A particularly interesting candidate is the pair density wave order, as an 
inherent mixture of density-wave and superconducting correlations {SI 5, S16, SIT). 

2.3 Simple nematic pseudogap order 

Nematic order, which breaks rotational symmetry but not translational symmetry of the lattice 
can also potentially explain our key observations in the superconducting state. We have stud- 
ied a simple form of a nematic distortion (Pomeranchuk type) of the bare Fermi surface that 
introduces a difference in t for the x and y directions, ty = at^ but not in t', t", t'" (eq is ad- 
justed to maintain the Luttinger's volume), a is chosen so that the experimental kc falls on 
the distorted Fermi surface and the superconducting gap opens along the distorted Fermi sur- 
face (inset of Fig. S|8t)). The tight-binding parameters are chosen to be tx,t' ,t'\t"\ eq & a = 
0.22, -0.034315, 0.035977, -0.0071637, -0.24127eV & 1.0989. 

Two types of orthogonal domains exist in a tetragonal crystal. ARPES is expected to detect a 
superposition of signals from both domains which gives rise to an apparent two-band band structure 
as shown in Fig. S|8K- The onset of superconductivity produces a gap for each band (Fig. 5(8^). 
Consistent with the experiment, two band features, one dispersionless and the other dispersive 
showing back-bending, are seen in the antinodal region (Fig. S(8]C). But two bands can still be 
seen away from this region (Fig. S[8j3), in contrast with what is seen in experiment. However, one 
should note that an energy-dependent broadening could in principle make two bands appear to be 
one when they are close enough. Where the two bands are farthest apart in the antinodal region, 
an increased coherence of the low-energy band at low temperatures could allow two features to be 
separately resolved. 
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2.4 Pitfalls 

Although these mean-field results present a favorable zeroth-order consistency with the experi- 
ment, it is important to point out that they also fail to capture other important aspects of the data. 
Most obviously, in mean-field theory the quasiparticles are exact eigenstates which produce sharp 
spectral features; this contrasts with the broad features seen in experiment, e.g., cf. Fig. S|7J-J 
near the M point. We have tried to correct this by including a phenomeno logical broadening (de- 
scribed below) but this still does not adequately reproduce the data (see the caption of Fig. $7]). 
These discrepancies probably reflect the effects of finite correlation length {S3) and strong electron 
correlations. 

2.5 Quasiparticle lifetime broadening introduced to the simulations 

All simulation results, except those in Figs. S|6]& S|7K, are shown with the renormalized band 
dispersion e'{k) broadened by a Fermi-liquid-like energy-dependent linewidth V = O.lco'^ + 0.005 
(eV) and its intensity proportional to Fig. S|3C assumes a phenomenological (marginal- 

Fermi-liquid-like) spectral function with a linewidth, F = ^ {au^ + {f3Ty (eV), where a =0.715 
and (3 =2.51 . This model was demonstrated to provide a reasonable global fit to the experimental 
antinodal spectra of heavily-overdoped Bi2212 {Tc= 65 K) taken at various photon energies (16.6 
- 32 eV) within energy (-200 - 50 meV) and temperature (10 - 128 K) ranges comparable to our 
study {S18). 

2.6 Effects of finite experimental resolutions and quasiparticle lifetime 

The ARPES simulation in Fig. [6]is based on a phenomenological spectral function in the supercon- 
ducting state (S19). We use the same global tight-binding bare band structure and assume d-wave 
superconductivity (A = 35 meV). The experimental energy and momentum resolutions and a re- 
alistic quasiparticle linewidth (a =0.715 and (3 =2.57) are incorporated. We find that these factors 
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have a negligible effect on the ka — kp misalignment. 

3 Additional discussion 

3.1 On the nature of the broken-symmetry pseudogap state 

An outstanding task after presenting our data in the main text is to identify the nature of the broken- 
symmetry state, and to relate it to the apparently similar electronic changes that occur below T* in 
other cuprates. There are two real-world considerations which complicate this program: 1) To the 
extent that the pseudogap state involves broken spatial symmetries, the role of quenched disorder 
is typically severe. It is expected to round the transition and to limit the growth of the pertinent 
correlations below T* to a limiting length scale which diverges only as the strength of the disorder 
tends to zero. This means that macroscopic measures of broken symmetry, especially where (as in 
Bi2201) the dopant atom distribution produces an inescapable degree of randomness, are expected 
to be less precise than they are in idealized models. 2) In comparing different materials, even 
if the same "universal" physics underlies the pseudogap, there can be detailed material-specific 
differences, even in the precise patterns of broken symmetry involved, making cross-material com- 
parisons more subtle than they would be in idealized models. 

While nominally the onset of PKE at T* suggests that the pseudogap phase uniformly breaks 
time-reversal symmetry, the small size of the effect and the lack of any directly detectable fer- 
romagnetism below T* make this conclusion uncertain, as previously discussed in the case of 
YBCO (SIO). A recent zero-field nuclear magnetic resonance (NMR) study detected no static 
magnetism on the Cu sites in superconducting Bi2201 (S20), implying that either the magnetism 
has a very novel character (S21, S22) or that there is no magnetism. In either case, the similarity of 
the PKE onset at T* in Bi2201 and YBCO supports the proposition that the pseudogap order has a 
common character in a broad class of cuprate superconductors. 

The large magnitude of the antinodal gap (~ S.S/csT*, measured relative to Ep) at tempera- 



11 



tures well below T* and the fact that it removes a substantial portion of the Fermi surface, implies 
that the spectral changes detected in ARPES reflect the primary pseudogap order. This is consis- 
tent with what has been inferred from the temperature dependence of the Knight shift measured 
by NMR (S20) of bulk Bi2201 in high enough magnetic fields to suppress superconductivity: The 
density of states (DOS) at Ep begins to drop sharply at T*, exhibiting a similar temperature depen- 
dence as those shown in Fig. 3. However, also consistent with the ARPES, a substantial fraction 
of the DOS remains ungapped down to K. 

3.2 Implications with various reported candidates for the pseudogap order 

Here we briefly compare the present results with studies of putative pseudogap order in other 
cuprate superconductors. Strong evidence has been found of nematic order in YBCO (Sll, S23) 
and Bi2212 (S24, S25). Evidence of unconventional translation-symmetry-preserving antiferro- 
magnetism in YBCO (S26) and HgBa2Cu04+5 (S27) has also been reported. These experiments 
suggest the pseudogap may primarily involve a subtle form of discrete symmetry-breaking with 
the wave vector Q = 0. While such order could produce significant shifts in the electron dis- 
persion, without breaking translational symmetry (or particle conservation), however, it is hard to 
see, at least at mean-field level, how it would lead to a gap in the spectrum. There is substantial 
evidence of unidirectional spin and charge density wave order (52S) in various cuprates, which 
spontaneously break lattice translational symmetry. Such order could readily lead to the opening 
of a gap on the nested portions of the Fermi surface. However, direct evidence of the existence 
of such a state has been found by neutron and X-ray scattering experiments only at temperatures 
below 50~ 60 K (well below T*) in particularly "stripe-friendly" materials. Unless the order 
parameter is somehow anomalously difficult to observe in diffraction (S21), it is hard to believe 
that similar experiments would not have already seen new Bragg peaks in the pseudogap regime 
of many cuprates. 
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3.3 On the ordering vector of a putative density- wave pseudogap order 

Ideally, by analyzing the momentum dependence of the spectral changes in ARPES spectra in- 
duced by the pseudogap, one might be able to back out information concerning the nature of the 
pseudogap order, such as the ordering vector, Q, of a putative density-wave state. In mean-field 
theory, starting from a bare band structure, e^, as the temperature T is reduced below T*, the 
pseudogap-induced spectral shifts would be expected to grow in proportion to (T* — T)^ with 
/3 = 1/2 only at special points in k space for which = e^_(^. In all other cases such shifts would 
be proportional to (T* — Tj^^/je-^ — From Fig. 3, to the extent that the spectral shift at kp 

and the change of Kerr rotation near T* can be associated with a critical exponent, they are both 
proportional to (T* — T)^^ with (3 ^ 1/2. The former is thus consistent with the mean-field expec- 
tation for a poor nesting of the antinodal portion of Fermi surface, which is in turn supported by 
the observed kc — kp misalignment. A nematic order with Q — expects a scaling of {T* — T)^, 
which seems not quite consistent with this particular aspect of our data. 

Although 0) can in principle be determined by studying the critical exponent as a function 
of position along the Fermi surface, there is every reason to think that order parameter fluctuations 
will significantly broaden both the temperature and momentum dependences of the spectroscopic 
signatures of the transition (even when the transition itself characterized by transport measurements 
remains sharp). Moreover, although the bare electronic structure is simple, its strange-metal nature 
precludes a safe description in terms of non-interacting quasiparticles with dispersion e^. 
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Figure S 1 : ARPES temperature dependence across T* A-E, Selected EDCs at five representa- 
tive temperatures along Cut CI shown in Fig. 1. The single EDC maximum is tracked to give the 
band dispersion summarized in the right inset of Fig. 3. F, Summary for the temperature-dependent 
binding energy position of the EDC maximum at M, kpi and kF2- Error bars are estimated based 
on the sharpness of features. 
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Figure S 2: Determination of the EDC shoulder feature by spectral subtraction. Subtraction 
was made for each EDC (e.g., the EDC at M at 10 K of Cut CI as shown in A) by assuming 
a linear background connecting intensity at Eji=-50 meV and Ej2=-10 meV. Intensity was then 
normalized by the linear intensity (I30) at Ej=-30 meV. The subtracted EDCs show a clear intensity 
maximum from which we define the position of the original EDC shoulder. This method yields a 
similarly non-dispersive peak at a similar energy position as the spectral division using the 40 K 
data (e.g., as shown in B for the same data set as in Fig. AC). It also provides a convenient way 
to quantify the temperature- and momentum-dependent existence of the EDC shoulder feature. C 
shows the subtracted EDCs (A) at various temperatures and E plots the temperature dependence of 
the maximum intensity (red) in comparison to a different subtraction using Eji=-60 meV and Ej2=0 
meV (blue). D & F show the results for Pl-PlO (in Fig. 2W). We defined phenomenologically 
the threshold for existence of the EDC shoulder feature at the half-drop positions of the maximum 
intensity w.r.t. the background (solid/empty symbols for existence/non-existence in E & F). This 
consistently supports the eyeball estimate that the EDC shoulder feature loses its clear definition 
above and beyond the antinodal region (along F-M). Note that its existence perpendicular to F- 
M could in principle extend wider than where is currently identified by green circles (Fig. 20-R), 
due to its potential mixing with the EDC maximum feature near back-bending. Also note that the 
non-vanishing variation of the maximum intensity above Tc in E is consistent with the idea that 
superconducting fluctations persist above Tc but still well below T*, in nearly OP Bi2201 (529) 
and other cuprates {S30, S31). 
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Figure S 3: Temperature dependence of time-resolved reflectivity. Results are shown at differ- 
ent laser powers, A 15 mW, BIO mW, and C 4 mW. The grey and pink shaded regions highlight 
the pseudogap and superconducting responses respectively. Note that the 15 mW data was taken 
over a shorter delay range than the lower power data. 
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Figure S 4: Gap function along the underlying Fermi surface at 10 K. A, Binding energy 
positions of various spectral features in the EDCs at kp plotted as a function of the Fermi surface 
angle. Red dashed curve is the gap function reported for OP Pb-Bi2201 (S32). Black dashed curve 
is a guide to the eye for a simple d-wave gap function. Our observation of multiple EDC features 
of comparable spectral weight in the antinodal region highlights the essential physics that eluded 
previous revelation. B, Selected EDCs atkp (kpi) as indicated by black arrows (dots) in (the inset 
of) A. The maximum feature is defined for EDCs with a single component seen away from the 
antinodal region. Both the shoulder at low energy and the hump at high energy are defined for 
EDCs at A;^;^ in the antinodal region. The (antinodal) EDCs at the neighboring ka (kai) are shown 
in dashed curves, which have a single component, reminiscent of those used to define the gap 
function in the antinodal region in many previous studies. 
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Figure S 5: Global tight- binding fit. A, Global tight-binding fit (red curves) to the EDC dis- 
persions at 172 K (black circles) whose momentum loci are indicated by blue lines in B. For 
clarity, adjacent cut has a 50 meV relative vertical offset. Because of the global constraint, the fit 
gives a slightly larger (~ 7 meV) band bottom at M point than the experiment. B shows that the 
Fermi surface given by the global fit coincides with the experimental one reproduced from Fig. 1. 
Please note that the experimental band dispersions are obtained (throughout this work) based on the 
EDC analysis, which is the traditional way of extracting dispersion relations from photoemission; 
it also suffers much less from the momentum-dependent matrix elements than the momentum- 
distribution-curve analysis, such that a parabola-like dispersion, rather than the "waterfall" type, 
remains to be seen close to the T point. 
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Figure S 6: ARPES simulation on the kc — kp misalignment issue. Even in a Bardeen-Cooper- 
Schrieffer (BCS) superconductor, 1) a finite shift in kp can result from a finite shift in the chemical 
potential assuming a particle number conservation in the gapped state. We find that the chemical 
potential increases by 5 meV in the superconducting state (A = 35 meV), which results in an 
effective bare band at low temperature (dashed black curve) with an increased kp relative to the 
bare band in the ungapped state (solid black curve). This increase is too small to be able to account 
for the kc ~ kp misalignment experimentally observed; 2) a finite kc — kp misalignment can in 
principle result from a finite momentum dependence of the superconducting order parameter along 
the cut. Again this effect alone is negligible in our case, as shown by the simulated false-color 
plot for the spectral function: the EDC peak dispersion (red dotted curve) shows a back-bending 
with kc very close to kp (of the solid black curve); 3) this simulation with realistic experiment 
conditions (see the Simulations part) also rules out a trivial interpretation of the observed ko — kp 
misalignment due to finite energy, momentum resolutions and quasiparticle lifetime. 
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Figure S 7 (following page): Simulations for the superconductivity-density- wave coexistence. 

A & E, Summary of experimental dispersions at 10 K and 172 K along Cut CI (reproduction of 
Fig. 20) and P1-P16 (Fig. 2V-W), respectively. B & F, C & G, D & H, Renormalized band disper- 
sions by simulations assuming different long-range orders: rf-wave superconductivity (SC, A = 35 
meV), its coexistence with bond-direction qi = (O.lSvr, 0) & g2 = (0, O.lSvr) checkerboard density 
wave (DW, Vi = 20 meV) or with bond-diagonal qaf = (tt, tt) density wave (V2 = 35 meV), re- 
spectively, perpendicular to or along F-M. Dashed curves are bare band dispersions resulting from 
a global tight-binding fit to the experimental EDC peak dispersions at 172 K (Fig. [S]). Renormal- 
ized band dispersions are moderately broadened for visualization purpose and the intensity reflects 
intrinsic spectral weight. Note that the visual separation between dominant states at low and high 
energies is not a real hybridization gap but due to suppression of intrinsic spectral weight of dis- 
persions within (C-D). Insets of E-H: kc (arrows in A-D) shown in half of a quadrant, from the 
experiment and simulations. Apparent discontinuities in the insets of G-H are due to changes of the 
physical character of back-bending. Error bars are smaller than the symbol size. I-K, Comparison 
between experimental and simulated EDCs at M point at 10 K and 172 K. Experimental intensity 
is normalized by the incident photon flux. Simulated EDCs in J are obtained from G, which are 
broadened phenomenologically in K (see the Simulations part). Different arrows are eyeguides 
for the EDC linewidth and energy shift of the EDC centroid. In K, multiple features below Ep at 
low temperature merge into a single feature, with its centroid energy and overall linewidth barely 
changed from high temperature, still very different from the experiment (I). L-N, Cartoons for the 
evolving band structure in different phases. Fermi surface is shown in black at Ep. Color scale is 
proportional to the intrinsic spectral weight in M & N. A (tt, tt) density wave order (V2 = 50 meV) 
and its coexistence with rf-wave superconductivity (A = 35) are assumed in M & N, respectively, 
only for illustrative purpose. Note that the values of parameters in all simulations are chosen only 
for a qualitative agreement with the experiment. 
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Figure S 8: Simulations for the superconductivity-nematic-order coexistence. A (B), Cartoon 
for the band structure in the nematic phase without (with) coexisting d-wave superconductivity 
(A = 35 meV). Two sets of bands from orthogonal domains are superimposed. C (D), Renormal- 
ized band dispersions perpendicular to (along) F-M in the coexisting state. Black (green) dashed 
curves are the bare (distorted) band dispersions. Renormalized band dispersions are moderately 
broadened for visualization purpose. Back-bending is pointed out by arrow. Insets: kc from the 
experiment that falls on the distorted Fermi surface in green. Color scale is proportional to the 
intrinsic spectral weight in B-C. Note that the values of parameters in the simulation are chosen 
only for a qualitative agreement with the experiment. 
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